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Compact binary coalescences, such as binary neutron stars or black holes, are among the most 
promising candidate sources for the current and future terrestrial gravitational-wave detectors. 
While such sources are best searched using matched template techniques and chirp template banks, 
integrating chirp signals from binaries over the entire Universe also leads to a gravitational-wave 
background (GWB). In this paper we systematically scan the parameter space for the binary coa- 
lescence GWB models, taking into account uncertainties in the star formation rate and in the delay 
time between the formation and coalescence of the binary, and we compare the computed GWB 
to the sensitivities of the second and third generation gravitational-wave detector networks. We 
find that second generation detectors are likely to detect the binary coalescence GWB, while the 
third generation detectors will probe most of the available parameter space. The binary coalescence 
GWB will, in fact, be a foreground for the third-generation detectors, potentially masking the GWB 
background due to cosmological sources. Accessing the cosmological GWB with third generation 
detectors will therefore require identification and subtraction of all inspiral signals from all binaries 
in the detectors' frequency band. 

PACS numbers: 



1. INTRODUCTION 

The ground-based gravitational-wave detectors are 
rapidly increasing their sensitivities. The first genera- 
tion detectors LIGO [l|, Q and Virgo 0, |3| have reached 
their design sensitivities and collected excellent data over 
several years of exposure. The second generation detec- 
tors, Advanced LIGO MM, Advanced Virgo 0, GEO- 
HF d, and LCGT d, [iB] are currently being built and 
commissioned. With 10 times better strain sensitivity, 
these detectors are expected to yield first direct detec- 
tions of gravitational- wave signals, and their first data is 
expected as early as 2014. Furthermore, there are already 
efforts under way to design the third-generation gravita- 
tional wave detectors, with another factor of 10 improve- 
ment in sensitivity. This includes the Einstein Telescope 
1^ . for which the design study was recently 
These detectors are expected to open a new 

providing new ob- 



have been performed using LIGO and Virgo data [15Hl 



It has also been argued that adding the gravitational- 
wave signals from all binaries in the Universe will produce 
a gravitational-wave background (GWB) - for example, 
see [IMQl for the most recent studies in the context 
of terrestrial detectors. The LIGO and Virgo collabora- 
tions have developed techniques for searching for GWB 
by cross-correlating data from pairs of gravitational wave 
25[- Such searche s have also been performed 



detectors 

using LIGO and Virgo data j26l - l28l |. and have produced 
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project 
completed. 

era in astronomy and astrophysics, 
servations of various events and objects in the Universe, 
complementary to the standard electromagnetic observa- 
tions. 

Among the many sources of gravitational waves, the 
coalescences of binary systems, such as binary neutron 
stars (BNS), binary black holes (BBH), or a black hole 
and a neutron star (BHNS) stand out as the most likely 
candidates for first detections. These systems gener- 
ate well understood "chirp" gravitational- wave signals, 
which have been computed using post-Newtonian ap- 
proximation [l3| or numerical relativity simulations [l4 | . 
One can then search for the chirp signals using matched 
template techniques - indeed a number of such searches 



competitive upper limits on the energy density carried 
by gravitational waves. 



The goal of this paper is to perform a detailed study 
of the accessibility of the GWB produced by the bi- 
nary coalescences to the second and third generation 
gravitational- wave detectors. Our study follows the work 
of Regimbau and Mandic [l3| , and includes detailed scans 
of the parameter space in these models, as well as possible 
effects due to the uncertainty in the star formation rate 
and in time-delays associated with the formation of the 
binaries. We will show that this background is likely to 
be observed by the network of second generation detec- 
tors, and that the third-generation detectors will likely 
be able to explore most of the parameter space for these 
models. In Section 2 we summarize the calculation of the 
energy density for these models. In Section 3 we present 
results of our systematic study, and we include conclud- 
ing remarks in Section 4. 



2. CALCULATION OF THE ENERGY 
SPECTRUM 

The energy spectrum of gravitational waves is usually 
described by the dimcnsionless parameter: 



/ dpGwif) 
Pc df 



(1) 



where / is frequency, dpcw/df is the energy density in 
the frequency range [/, f+df] and pc is the critical energy 
density needed to close the Universe: 
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where Hq and G are the Hubble parameter and Newton's 
constant respectively and c is the speed of light. 

The energy spectrum for the case of binary coales- 
cences can be written as follows (see for instance [3; SH- 
23[, we will follow ^19|] in our approach): 



f^Gw(/) = —fF{f) 

PcC 



(3) 



where the integrated flux (per unit frequency) is defined 
as: 
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where Rz{z) is the rate of gravitational- wave sources per 
interval of redshift z as observed in the detector (Earth) 
frame, dL{z) — (l + z)r(z) is the luminosity distance, 
r(z) is the proper distance, and is the gravitational 

spectral energy emitted by a single source and observed 
in detector frame. The rate in Eq. 2] is given by: 



Rz{z) = XRviz) 



dVjz) 
dz 



(5) 



where A is the mass fraction converted into progenitors 
(discussed in more detail below), Rv{z) is the observed 
rate of binary coalescences (in units of mass per unit 
comoving volume per time), and 



dV{z) _ Anc r2(z) 



dz Hq E{D,m,^a,z) 



(6) 



with E{Q,M,^A, z) = ^ynMi^ + z)^ + capturing the 
dependence of the comoving volume on redshift. We use 
the standard ACDM cosmology, with flu = 0.3, JIa — 
0.7, and Hubble parameter Hq — 70 km s^^ Mpc^^. The 
rate Rv{z) is dependent on both the star formation rate 
and on the time-delay td between the formation of the 
binary system and the actual coalescence, and can be 
written in the following form: 



Rv{z) = 



1 



1 + Zi 



R.{t,{z) ~ td)P{td)dtd (7) 



where i?* is the star formation rate (discussed further 
below), tc{z) is the cosmic time corresponding to redshift 
z, and Zf \s the redshift at the formation time tc{z) — td- 
The factor {1 + Zf) in the denominator corrects for the 
time dilation due to the cosmic expansion and converts 
the rate from the source frame into the detector frame. 

Population synthesis [2^, [s^] suggests that the prob- 
ability distribution for the delay time is well described 
by P{td) ~ for td > imin, where imin is the minimum 
delay time for a massive binary to evolve until coales- 
cence. While the currently preferred parameter values 
are a = -1 and imin = 20 Myr for BNS and 100 Myr for 
BBH, other values cannot be excluded. Following [3l|. 
we will investigate the following ranges for these param- 
eters: a = -0.5, -1,-1.5, imin = 20, 100 Myr for BNS, 
and imin = 100, 500 Myr for BBH. We will also examine 
the case where time delay is ignored, td = 0. The star 
formation rate i?* has also been investigated by several 



authors 32 



2i,aa. We win 



-36j, for a recent review see 
investigate dependence of our results on the choice of the 
star-formation rate. 

The scaling factor A (in units of Mq^) is a parameter 
which includes three different effects: for BNS, these are 
the mass fraction of neutron star progenitors, the fraction 
of massive binaries formed among all stars, and the frac- 
tion of binaries that remain bounded after the second su- 
pernova event (and similarly for the BBH and BHNS). All 
of these factors are associated with significant uncertain- 
ties, which is why we will treat A as a free parameter of 
the model in our study. We note, however, that Ai?y(0) 
represents the local (present) rate of binary coalescences. 
These local rates have been a subject of multiple stud- 
ies, as they directly impact the number of individual bi- 
nary coalescences that could be detected by the second- 
generation gravitational-wave detectors. A recent study 
by the LIGO Scientific Collaboration has produced pes- 
simistic, realistic, optimistic, and maximal possible esti- 
mates for these rates fsij , based on the observed galactic 
binary pulsars and on the population-synthesis models. 
We will compare the results of our study to these rates 
estimates. 

The final factor appearing in the Equation [4] is 
dEQ^^/df, the gravitational spectral energy from a sin- 
gle source. For the BNS and BHNS models, we will only 
include the inspiral part of the gravitational- wave signal. 
In the quadrupolar approximation, and assuming circu- 
lar orbit, the observed spectral GW energy, averaged over 
orientation, from a binary system at redshift z is given 
by a rather simple form: 



df 



2/3 



{M^'f^f-^/^ 



(8) 



where = (1 -I- z)Mc is the observed redshifted chirp 
mass of the binary system and Ale is the physical mass. 
We will assume the following ranges for Mc- 1-2.5 Mq 
for BNS, 2.5-10 Mq for the BHNS, and 2.5-20 Mq for 
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the BBH models. These mass ranges are based on [30j, 
allowing for uncertainties in possible neutron star and 
black hole masses. For the BBH case, however, we will 
follow (isl . [2^ and use the more complex functional form 
derived by [24| . which includes the inspiral, merger, and 
ringdown contributions to the gravitational-wave signal. 

The upper limit on the integral range in Eq. |4] depends 
on both the emission frequency range, /min — /max, hi 
the source frame, and on the maximum redshift Zmax 
considered for the star formation history calculation: 
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Combining the expressions above, we obtain for the 
density parameter: 



^gw (/) — 
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(10) 

Unless noted otherwise, we set Zjnt = in our calculation. 

We emphasize, however, that the GWB computed here 
is not necessarily continuous in time, as already noted in 
[H H m ill. To illustrate this we compute a duty 
cycle parameter, defined as 



Rz (z) dz 



(11) 



where dT(z)/df represents the time a binary spends in 
the frequency band [/, / + df] after properly accounting 
for redshift: 



dr _ be' 
'df ~ 967r8/3G5/3 



(12) 



The quantity dA/df then captures the number of bina- 
ries generating gravitational-wave signals in a 1 Hz bin 
as observed by a detector on Earth (it can be compared 
to the overlap function of 2^, which is the number of 
sources present on average in a frequency bin A/ around 
the frequency /). This quantity is plotted in Figure [T] 
for three different types of binaries and for different val- 
ues of Zsup- Note that integrating the BNS, Zsup — 6 
curve over the Advanced LIGO frequency band (roughly 
10-200 Hz) yields A ^ 10 - in other words, in any 0.1 
sec long time-segment (corresponding to the lowest ob- 
servable frequency of 10 Hz) there will be on average 10 
binary neutron star systems emitting in the 10-200 Hz 
band. The duty cycle is somewhat lower for the BBH 
and BHNS cases. 

Figure [T] also shows that most of the contributing bi- 
naries reside at redshifts z > 0.1. The nearest bina- 
ries are expected to produce loud chirp-like signals that 
could be individually detected by the upcoming detec- 
tors. However, such loud transients are typically explic- 
itly excluded from the searches for GWB |25l428| . We 



have verified that the nearest (and loudest) binaries con- 
tribute little to ^Igw ■ Figure [2] shows the gravitational- 
wave spectrum ^Igwif) computed for the BNS case with 
Mc = I.22M0, A = 3 X IO^^Mq^ star formation rate 
from and P{td) ^ t^^ with i^in = 20 Myr. Ex- 
cluding the nearest binaries (e.g. those with redshifts 
z < 0.1) leads to a small (< 2x) reduction in the spec- 
trum amplitude. We have further verified this conclu- 
sion with explicit Monte Carlo simulations for the case 
of Advanced LIGO collocated detector pair. Note that 
a similar Monte Carlo simulation was performed in the 
context of the Einstein Telescope |4fl] . 
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FIG. 1: Number of binaries per 1 Hz frequency bin. For the 
BNS cases we assume each star to have mass of 1.4 Mq, local 
rate of A = 1 Mpc"^ Myr"\ and tmin = 20 Myr {P{td) ~ 
1/t). For the BHNS case we assume masses of 1.4 Mq and 10 
Mp,, the local rate A = 0.03 Mpc"^ Myr~\ and t^in = 100 



^0, 

Myr. For the BBH case we assume masses of 10 Mq, 
W„ = 100 Myr. 



and 



3. RESULTS 

We performed a systematic study of the GWB due to 
binary coalescences, described in Section 2. In particular, 
we performed a scan of the parameter space spanned by 
the parameters A and Mc for each of the BNS, BBH, and 
BHNS cases. We compare the model predictions for each 
point in the parameter space with sensitivities of current 
and future gravitational wave detectors, as well as with 
the estimates of the local coalescence rates presented in 
[sst . We also investigate the importance of the choice of 
star formation rate and of the choice of the probability 
distribution for the time delay between the formation of 
a binary and its coalescence. 

To start with, we investigate the redshift dependence 
of the integrand in Equation [10] for different choices of 
the star formation rate. More specifically, the integrand 
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FIG. 2; Gravitational wave spectrum i^g^if) computed for 
the BNS case with Mc = 1.22Mq, A = 3 x 1O"'^M0\ star 
formation rate from [s^l, and P{td) ~ with tmin = 20 Myr. 
The effect of removing the nearest binaries (with redshifts 
z < 0.1 or z < 0.2) is very small. 



can be written as: 
/(z) = 
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E{nM,nA,z){i + zy/^- 

Figure |3] shows I{z) for five different choices of the star 
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FIG. 3: I(z) for different star formation rates, assuming a — 
-1 and tmin = 20 Myr. 

formation rate. Since there are non-negligible differences 
between these five estimates of the star formation rate, 
we will present the results for the two extreme cases, 
namely the Hopkins & Beacom [sij and Nagamine et al 
(35| . In particular, we scan the parameter space in the 
A — Mc plane: for each point in this parameter space, we 
compute ^gw(f) and we compare it to the most recent 
upper limit from LIGO [i^ and to the projected sensi- 
tivities for the second-generation (Advanced LIGO) and 
the third-generation (Einstein Telescope) gravitational- 
wave detectors. Figure H] shows the results of the scan of 



the A — Mc plane for the two estimates of the star for- 
mation rate and for the three binary coalescence cases: 
BNS, BBH, and BHNS. For each of these cases, we ob- 
serve that the latest GWB upper limit obtained using the 
LIGO data |26i] excludes the largest values of A, which are 
larg er than the maximal expected local coalescence rates 
|38| . However, the values of A corresponding to the opti- 
mistic and realistic (in the case of BNS and BHNS) local 
coalescence rates will be accessible to the second gen- 
eration network of gravitational wave detectors (assum- 
^tandard Advanced LIGO expected strain sensitivity 
for two collocated detectors with one year of expo- 
sure). The values of A corresponding to the pessimistic 
local coalescence rates will be accessible to the third- 
generation gravitational-wave detector network (assum- 
ing ET-D strain sensitivity curve [ill, for two collo- 
cated detectors and one year of exposure). In fact, the 
binary coalescence GWB will be a foreground masking 
the GWB background due to early-Universe sources (in- 



flationary models [4lJ, |42|, or phase transitions models 
[43! ) which may be one of the targets of the third gener- 
ation detectors. Hence, to detect the cosmological GWB 
it will be necessary to identify and subtract all of the in- 
spiral signals from all binaries in the frequency band of 
the third generation detectors. This is a daunting task, 
but appears to be plausible as demonstrated in |44] for 
the framework of the Big Bang Observer satellite-based 
detector . These conclusions hold for all of the three 
cases (BNS, BBH, and BHNS) and are rather weakly de- 
pendent on the choice of the star formation rate (we have 
verified that using the remaining three estimates of the 
star formation rate [33|, [SJ, [3^] yields results roughly be- 



32| 



tween those shown here for the Hopkins & Beacom 
and Nagamine et al fiBJ). 

We also investigated the effect of different choices of 
the probability distribution P{td)- As noted above, the 
population synthesis suggest P{td) ^ for td > tmin, 
where tmin is the minimum delay time for a massive bi- 
nary to evolve until coalescence. Since there is some un- 
certainty in the parameters a and tmin, we will probe the 
range of values of these parameters discussed in the lit- 
erature. In particular, we examine a = —0.5, —1, —1.5, 
= 20, 100 Myr for BNS and tmin = 100, 500 Myr 
for BBH, as well as the case when there is no time de- 
lay between the formation and coalescence of the binary. 
Figure [5] shows the variation in the contours for the sec- 
ond and third generation detectors, for BNS and BBH 
models, for several of the P{td) parametrizations. We 
observe that the contours are rather insensitive to P{td), 
varying by at most a factor of 2 in the A parameter. 
Hence, the choice of P{td) does not qualitatively affect 
the conclusions of this study. 

Finally, we note that for the case of BBH a similar 
study was performed in - they computed Advanced 
LIGO and ET-D contours in the A — Mc plane that were 
substantially higher in A (~ 20 x for Advanced LIGO). 
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FIG. 4: Accessibility of binary coalescence GWB to current and future gravitational wave detectors. The two columns corre- 
spond to two estimates of the star formation rate: Hopkins & Beacom [33] (left) and Nagamine et al [s^ (right). The three 
rows correspond to BNS, BBH, and BHNS respectively, top to bottom. For each plot we show the A — Mc plane: the region 
of the parameter space excluded by the S5 LIGO result [2a| , and the expected sensitivities of the Advanced LIGO collocated 
detector pair (assuming 1 year of exposure B, [Hi), and of the Einstein Telescope (assuming two collocated detectors with ET-D 
sensitivity and one year of exposure [Tlji UM)- These regions are to be compared with the expected local coalescence rates 
shown as horizontal dashed lines: top-to-bottom they correspond to maximal, optimistic, realistic, and pessimistic estimates 
presented in [3^ . 
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These differences largely come from the different assump- 
tions in detector sensitivity and exposure. For Advanced 
LIGO they assumed a non-standard detector strain sen- 
sitivity (this is the dominant cause of discrepancy), non- 
collocated detectors, and 3 years of exposure - we assume 
the standard strain sensitivity, collocated detectors, and 
one year of exposure. For ET-D, they assumed the tri- 
angular detector configuration (leading to the factor of 
3/8 in overlap reduction), while we assumed L-shaped 
interferometers . 



4. CONCLUSIONS 

In this paper we computed the gravitational wave back- 
ground due to coalescences of binary neutron stars, bi- 
nary black holes, and black hole - neutron star binaries, 
following the approach of ^IB] • While such computations 
have been done in the past, in this study we performed 
a systematic scan of the parameter space, taking into 
account the possible variations in the result due to the 
choice of the star formation rate, and due to the choice 
of the distribution P{td) of the delay time between the 
formation and coalescence of the binary. For each point 
in the parameter space, we compare the model prediction 
to the expected sensitivities of the second and third gen- 
eration gravitational-wave detector networks. We find 
that models corresponding to the optimistic and realistic 
(in the case of BNS and BHNS models) local coalescence 
rates will be accessible to the second generation detec- 
tor network. We also find that models corresponding to 
the pessimistic local coalescence rates will be accessible 
to the third-generation detector network. The binary 
coalescence GWB will, in fact, be a foreground for the 
third-generation detectors, and it will mask the GWB 
background due to early-Universe sources. Accessing the 
cosmological GWB with third generation detectors will 
therefore require identification and subtraction of all in- 
spiral signals from all binaries in the relevant frequency 
band. 
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